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I.  INTRODUCTION 


Apparently  no  method  is  readily  available  in  the  published  literature 
which  will  enable  quick  and  easy  evaluation  of  the  total  heat  transfer  to  a 
body  in  a hypersonic  flow  environment.  Such  information  is  very  useful, 
for  example,  in  evaluating  heat  sink  or  internal  coolant  requirements  for 
various  body  shapes  in  a ground  test  hypersonic  environment,  or  the  trans- 
piration-cooled nosetip  coolant  and  reentry  vehicle  heat  shield  requirements 
for  ballistic  hypersonic  reentry  vehicles. 

The  initial  work  in  this  area  is  that  of  Allen  and  Eggers  (Ref.  1). 

These  authors  give  the  form  of  an  equation  for  total  heat  load  to  a reentry 
body.  However,  direct  calculation  of  results  is  not  possible  since  the  con- 
vective heat  transfer  relationship  is  left  unspecified  in  the  final  equation. 
Brunner  (Refs.  2,  3)  has  extended  the  work  of  Allen  and  Eggers  by  incorpo- 
rating flat-plate-reference-enthalpy  convective  heat  transfer  expressions  for 
both  laminar  and  turbulent  flow.  For  a complete  reentry  trajectory, 

2 

Brunner's  work  allows  the  evaluation  of  the  total  local  heat  load  in  Btu/ft 
at  a given  station  on  a sphere-cone.  However,  to  obtain  the  total  heating  to 
the  entire  body  from  Brunner's  results,  an  integration  must  be  performed  of 
the  given  laminar  or  turbulent  heat  transfer  distribution  functions  over  the 
sphere,  or  sphere-cone  body,  of  interest.  Also,  for  sphere-cones,  Brunner's 
results  are  approximate  because  of  the  approximate  pressure  distribution 
used. 

^Allen,  H.  Julian  and  A.J.  Eggers,  Jr.,  "A  Study  of  the  Motion  and  Aero- 
dynamic Heating  of  Ballistic  Missiles  Entering  the  Earth's  Atmosphere  at 
High  Supersonic  Speeds,  " NACA  Report  1381,  1958. 

2 

Brunner,  M.  J.,  ''Analysis  of  the  Aerodynamic  Heating  for  a Re-entrant 
Space  Vehicle,  " J.  of  Heat  Transfer,  August  1959,  pp.  223-229. 

3 

Brunner,  M.  J.,  "Heat  Protection  System  Concept  Evaluation,  " Report 
PIR  U-91 90-73- 1 14,  General  Electric  Company,  Philadelphia,  Pennsylvania, 
9 January  1973. 


The  primary  purpose  of  this  work  has  been  to  provide  analytic  equations 
and  calculated  numerical  results  which,  when  taken  together,  allow  the  total 
heat  transfer  for  body  shapes  of  interest  to  be  quickly  and  accurately  calcu- 
lated. 

Total  heat  transfer  relationships  for  both  constant,  ground  test  type, 
environments  and  realistic  ballistic  trajectory  environments  are  derived  in 
Section  II.  In  the  later  sections,  numerically  calculated  results  for  the 
laminar  and  turbulent  heat  transfer  integrals  are  given  for  families  of 
sphere-cones  and  flat-face  cones.  When  the  numerical  results  are  combined 
with  the  analytic  equations,  total  heat  transfer  in  both  types  of  environments 
may  be  readily  calculated  for  a wide  range  of  sphere-cone  and  flat-face  cone 
geometries. 
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II.  INTEGRATED  HEAT  TRANSFER  EXPRESSIONS 


A.  LOCAL  HEAT  TRANSFER  RATES 

An  expression  is  needed  for  the  total  integrated  heat  transfer  in  hyper- 
sonic flow  environments  as  a function  of  appropriate  flowfield,  body  scale, 
and  body  geometry  parameters.  But  first,  an  appropriate  expression  is 
required  for  the  local  heat  transfer  coefficient.  In  Refs.  4 and  5,  the  boundary 
layer  integral  momentum  equation  is  solved,  and  Reynolds  analogy  and  a 
compressibility  transformation  are  applied  to  obtain  the  following  expressions 
for  the  local  heat  transfer  coefficient: 


These  expressions  have  been  found  to  agree  very  closely  with  experimental 

data  at  high  supersonic  Mach  numbers  (Ref.  6);  limited  comparison  with 

% 

experimental  data  for  a hypersonic  Mach  number  also  shows  good  agreement. 

4 

Vaglio- Laurin,  R.,  "Laminar  Heat  Transfer  on  Three-Dimensional  Blunt 
Nosed  Bodies  in  Hypersonic  Flow,  " ARS  Journal,  29(2),  February  1959, 
pp.  123-129. 

5 

Vaglio- Laurin,  R.,  "Turbulent  Heat  Transfer  on  Blunt  Nosed  Bodies  in  Two- 
Dimensional  and  General  Three-Dimensional  Hypersonic  Flow,  " J.  of  the 
Aero/Space  Sciences,  27(1),  January  I960,  pp.  27-36. 

L 

» Widhopf,  G.  F.  and  R.  Hall,  "Transitional  and  Turbulent  Heat  Transfer 

Measurements  on  a Yawed  Blunt  Conical  Nose  Tip,  " ALA  A Journal,  10(10), 
October  1972,  pp.  1318-1325. 

. # 

Crowell,  P.  G.,  private  communication,  January  1973. 


-7- 


If  we  define  the  following  nondimensional  quantities  denoted  by  the  bars 


P 


e 


(2) 


then  the  local  heat  transfer  coefficient  expressions  become 


C.r  nP 
1 e 


'H 


/_  n+1  _ 
p u p r ds 
e e e 


n/(n+l) 


1/2 

'B 


n/  (n+1) 


(3) 


where  for  laminar  flow  Cj  = 0.  332,  n * 1 and  for  turbulent  flow  = 0.0288, 
n = 1/4. 

We  now  distinguish  between  a constant  freestream  environment  problem 
for  which  integration  over  the  surface  are-  gives  the  total  heat  transfer  rate 
in  Btu/sec  and  a trajectory  environment  problem  for  which  integration  over 
the  surface  area  and  for  all  altitudes  gives  the  total  heat  transfer  in  Btu's. 
This  is  done  below. 

B.  CONSTANT  FREESTREAM  ENVIRONMENT 

To  obtain  the  total  heat  transfer  rate,  we  integrate  the  local  heat  trans- 
fer rate  q over  the  surface  area  of  the  body.  Thus 


dQ 

dt 


(4) 


The  total  heat  transfer  rate  then  becomes  upon  substituting  from  Eq.  (1)  and 
nondimens  ionalizing 


dQ/dt 


OO  00  o\  Dt 


,(*4) 


■ ¥)l 


fl.%CHrd" 


(5) 
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T 


The  hypersonic  approximations  J 2Hq  «V  and  pj p^  = (y+l)/(y-l)  have  been 
made  in  obtaining  Eq.  (5).  When  Eq.  (3)  is  substituted  into  Eq.  (5) 

s 


dQ/dt 


2 C 


im\ 

i \ y- 1 / 


P V H [ReJ 

GO  oo  o \ B/  o 


n/ (n+1) 


/ 


— — 77  —n  + 1 
p \x  p r ds 
e e e 


/P  n P 
e e 


-i  n/  (n  + 1) 


" - 77  7n+1ds 

e 


(6) 


where 


Re  = 
o 


*,o(2Ho)1/2rB 

P„ 


(7) 


jjc 

As  noted  in  Ref.  7 (and  footnote  below  ),  the  integral  in  Eq.  (6)  is  of  the  form 
fil Q]'n/(n  + 1)dIQ  where 


— — t;  — n+1 
p u p r ds 
e e e 


(8) 


Thus,  Eq.  (6)  may  be  integrated  once  to  obtain 


dQ/dt 

p V H (nrt) 
oo  oo  o \ B / 


2<n+l>Cl(£i)  _ 1/(„+1) 

[R.  Q 

L oJ 


(9) 


Substituting  the  appropriate  values  for  and  n,  we  obtain 

(Laminar)  dQ.L/dt  - = 1.328  (ffi)  (i^ ) ^Re^72  (10a) 

P v H (nr-.) 

oo  oo  o\  B/ 


(Turbulent)  dQ  /dt  y—  - 0. 0720  (y±j)  (il ) Re‘1/5  d°b) 

p V H (ml)  'TnV  ' Q' 

oo  oo  o \ B/ 


7 

Aihara,  Y.,  "Optimum  Body  Geometries  of  Minimum  Heat  Transfer  at 

Hypersonic  Speeds, 7 * * *  11  ALA  A Journal,  Technical  Notes,  6^(11),  November 

1968,  pp.  2187-2188. 

* Crowell,  P.G.,  private  communication,  November  1976. 
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c. 


TRAJECTORY  ENVIRONMENT- - BASIC  EQUATIONS 


In  this  case  to  obtain  the  total  heat  transfer,  we  must  not  only  integrate 
the  local  heat  transfer  rate  over  the  surface  area  of  the  body,  but  must  also 
integrate  in  time  (or  altitude)  through  the  trajectory.  If  we  let  £/€-  total  heat 
transfer /unit  area,  then 


dJtf  • 
~dt  = q 


(ID 


The  development  here  closely  follows  that  given  by  Allen  and  Eggers  (Ref.  1) 
except  that  a key  simplifying  assumption  has  been  eliminated;  i.e.,  Allen 
and  Eggers  assumed  that  Re^  is  not  a function  of  altitude.  If  altitude  is 
denoted  by  y,  to  a good  approximation 

= -V  sine,.  (12) 

dt  oo  E 

also 

§’J f/Wds  (13) 

S 


After  Eqs . (11)  through  (13)  are  combined  and  integrated  circumferentially 


dQ 

dy 


-qrds 

V sin  6 
oo 


E 


(14) 


Proceeding  with  nondimensionalization,  approximations  and  simplifications 
exactly  as  in  the  constant  environment  case  above,  we  obtain 


dQ 

W - 


-',a»V®’rrB 


sin#. 


u C..r  ds 
e e H 
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(15) 


Substituting  for  CV  from  Eq.  (3),  we  find 


dQ 

dy 


.C.p  V2 

1 00  OP  BVy-1/ 

[Re 

L o E 


— — — — n+1  — 
p u U r ds 


■ s 

f — — 77  —n+1 

1 p u H r ds 

J 'e  e^e 

n/(n+l) 

L 0 

(16) 


When  we  integrate  in  s as  before 


dQ 

dy 


-(n+l)C1 


n ,.2  2 /y+i\ 

PooVoo7rrBlFT) 


sinVReo] 


n/(n+l) 


■r  1/ (n  + 1) 
Q 


(17) 


Before  we  can  integrate  in  y,  we  must  express  p^,  and  ReQ  as  functions 
of  y.  Notice  that  because  of  the  way  in  which  it  is  nondimensionalized,  Iq  is 
not  a function  of  y. 

To  a very  good  approximation,  the  viscosity  over  the  range  of  pressure 
and  temperature  of  interest  may  be  assumed  to  be  independent  of  pressure 
and  represented  by  a power  law  temperature  dependence.  Thus 


"o  - Vo  <18> 

When  we  use  ideal  gas  expressions  and  employ  the  hypersonic  approximation 
72^  = V , Eq.  (18)  becomes 


"o 


v2  0) 
oo 


Substituting  into  Eq.  (7),  we  may  write  ReQ  as 


(19) 


(20a) 





The  refore 


Ren/(n+1)  = 
o 


n / (n+1) 

P co 

y(2U-l)n/  (n+1) 
/ oo 


(20b) 


where 


n/  (n+1) 


(21) 


From  Allen  and  Eggers  work 


P 


oo 


V 

00 


= p0.-** 


* VEe 


B e-Xy 

T 


(22a) 

(22b) 


where  the  ballistic  factor  B and  the  ballistic  coefficient  £ are  given  by 


B - , f° 

X0sin#E 


(23a) 


0 = 


m 


cda 


(23b) 


Substituting  Eqs.  (20b)  and  (22)  into  Eq.  (17),  we  obtain  finally 


f(2o>+l)n+2]B  -Xy 
dQ  _ _ -Xy/(n+l)  ” 2(n+l) 

37  c3 


(24) 


where 


C 


3 


(n+lJCj 

nrB(^r) 

C2 

sin0E 

l/(n+l)  [(2w+l)n+2]/(n+l)  l/(n+l) 

’°  VE  *Q 

(25) 
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D.  TRAJECTORY  ENVIRONMENT-- LAMINAR  FLOW 

When  we  substitute  for  n,  C^,  and  in  Eqs.  (24)  and  (25) 


dQ^  _ „L  -(A/2)y  - (2W+3)B/4  e'Xy 

dy  ' C3  e 


(26a) 


* -0.939 


r /y-1  \ufY+l\ 

^n\zyR/  \ y-  1/ 


B 

sin0r 


1/2  (2UH-3)/2  f "1 1/2 


P°  V, 


IL_ 

. Q 


(26b) 


When  the  limits  of  integration  are  reversed  in  altitude  y 


qL  = _CL  C e-(X/2)y  e-[(2W+3)B/4]  e' 


Define  a laminar  ballistic  factor  B^  as 


lL  = <2“±-3)  B 


Equation  (27)  can  then  be  integrated  directly  to  obtain 


MbM1 


y(i/2,  bl> 


where  y(l/2,  B ) is  the  incomplete  gamma  function. 

Upon  substitution  for  and  the  definition  of  the  incomplete  gamma 
function,  Eq.  (29)  may  be  rewritten  as 


QL  = 0.  93  9 


H 0°  v 
Xsin0E  E 


1/2  (2w+3)/2  _ .1/2 


(£)  -4^) 


-13- 


Finally,  substituting  for  and  B from  Eqs.  (28)  and  (23a),  we  then  have 


QL  = FL  erf  (/b17  ) 


FL=  1.  866x  10'6 


Ml 

(2w+3) 


'y+i \ I ' c /y-i\  5il  ; 

SUyRj  _ 3 / 2 „ (3+2(j)  /2 

Asin0  rB  V E 


(■T 


e * dt 


(31a) 


(31b) 


(31c) 


E.  TRAJECTORY  ENVIRONMENT- - TURBULENT  FLOW 

When  we  substitute  for  n,  C . , and  C7  in  Eqs.  (24)  and  (25) 


dQT  _T  -(4X/5)y  -[(2w+9)B/10]e'Ay 
8 C3  e 


(32a) 


= -0.0627 


: (2=1.)' 
n\zyR> 


cj  -i  1 /5 


.(fS) 


4/5  (2<J+9)/5r  -«4/5 

^ VE  5 


(32b) 


Again,  if  the  limits  of  integration  in  altitude  y are  reversed 


~t  „t  r 

Q = -C3  I e 


■(4X/5)y  -[  (2w+9)B/  10]e_Ay  . 


Define  a turbulent  ballistic  factor  B 


T (2w+9)r 
B * ~ Tq— B 
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Integrating  Eq.  (33),  we  obtain 


» 

I 

* 


I 


I 


I ! 


i 


C 

QT  = - A-nr  y(4/5,  BT)  (35) 

X[bt]4/5 

Substituting  for  , B1  and  B as  in  the  laminar  case,  we  obtain 

QT  = FT  y(4/5,  BT)  (36a) 


FT=4.  70x  10-6 


(Z±i) 1 

4/5 

fc  (y-1  H 

\y-i/ 

^l2?R/ 

(2xH-9) 

Xsindg 

1/5 


v (9+2<J)/5 
E 


(36b) 


T 

where  y( 4/5,  B ) is  the  incomplete  gamma  function  given  by 


y(4/5,  BT) 


•'A 


rt  ,4/5 
B } .5/4 

e dt 


(37) 


Integrations  entirely  analogous  to  those  performed  on  Eqs.  (27)  and 
(33),  for  laminar  and  turbulent  flow,  respectively,  were  carried  out  by 
Brunner  (Refs.  2,3).  Brunner's  starting  point,  however,  was  the  "refer- 
ence enthalpy"  heat  transfer  expressions  for  laminar  and  turbulent  flow 
rather  than  the  Vaglio- Laurin  expressions  used  herein.  Also,  Brunner's 

final  results  for  integrated  heat  transfer  through  a trajectory  were  given  in 

2 

terms  of  Jf,  i.e.,  the  time  integrated  heat  transfer  in  Btu/ft  at  a point  on 
the  body. 


* 


Proper  units  for  Eqs. 


(31a)  and  (36a)  are  discussed  on  page  35. 
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III.  CALCULATED  RESULTS--CONSTANT 
FREESTREAM  ENVIRONMENT 


If  the  free  stream  environment  is  not  changing  with  time,  then  the  total 
heat  transfer  rate  can  be  expressed  in  terms  of  the  freestream  flowfield 
variables  and  the  body  geometry  as  given  by  Eqs.  (10a)  and  (10b).  These 
expressions  are  shown  graphically  in  Fig.  1.  The  behavior  is  seen  to  be 
similar  to  classical  Nussett  number-Reynolds  number  behavior,  with  a 
0.  5 power  Reynolds  number  dependence  in  the  laminar  case  and  a 0.  2 power 
Reynolds  number  dependence  in  the  turbulent  case.  Notice  that  the  right- 
hand  sides  of  Eqs.  (10a)  and  (10b)  can  be  interpreted  as  representing  the 
fraction  of  the  intercepted  freestream  total  energy  flux  which  is  transferred 
as  heat  to  the  body. 

Figure  1 is  applicable  for  all  body  shapes.  However,  the  ordinate 
contains  the  heat  transfer  integral  Iq,  which  depends  upon  the  shape  contour 
of  the  body  and  is  given  by  Eq.  (8). 


The  evaluation  of  In  reduces  to  evaluating  the  local  pressure  p at 
each  position  on  the  body  and  then  integrating  Eq.  (8)  because  the  normalized 
product  P&Pe  can  be  written  in  terms  of  the  normalized  pressure  pg  = Pe/Pt  as 


P P 
e e 


P -M-  _ [(y-i)0Hl]/y 

e --  = p 

P P *e 
roo 


(38) 


Also,  assuming  an  isentropic  expansion,  the  local  velocity  ug  can  be  written 

in  terms  of  p as 
e 


u 

e 


- (r-i)/y 

* m 


(39) 


Thus,  Iq  is  a function  of  the  specific  heat  ratio  Y,  the  temperature  exponent 
in  the  viscosity  law  to,  and  body  geometry. 


The  inviscid  flow  pressure  distributions  needed  to  determine  the 
boundary  layer  edge  properties  as  given  by  Eqs.  (38)  and  (39)  were  obtained 
from  the  best  numerical  inviscid  flowfield  calculation  procedures  available 
(Refs.  8-10).  A schematic  diagram  of  the  inviscid  flowfield  calculation 
procedure  for  flat-face  cones  is  shown  in  Fig.  2.  For  sphere-cones,  the 
start  line  for  the  rotational  characteristics  computer  program  (Ref.  9 ) was 
obtained  using  the  Ames  blunt  body  computer  program  (Ref.  8). 


The  basic  body  geometry  considered  in  this  work  is  a cone.  Calcula- 
tions were  carried  out  for  cone  half-angles  ( 6 ) of  30,  20,  10,  5 and  0 deg 
for  both  hemispherical  and  flat-face  noses.  In  each  instance,  calculations 
were  made  by  numerically  integrating  Eq.  (8)  starting  at  the  stagnation  point. 
Results  of  the  numerical  inviscid  flowfield  calculations  for  these  geometries 
were  used  in  conjunction  with  Eqs.  (38)  and  (39)  to  evaluate  the  integrand 
along  the  body.  In  this  way,  Iq  and  Iq  for  each  cone  half-angle  and  a con- 
tinuous range  of  thickness  ratio,  T = Zr^/lj,  were  calculated.  The  bluntness 

ratio,  B.  R.  = and  the  thickness  ratio  are  related  by  Eqs.  (40a)  and 

In  r> 

(40b).  For  sphere-cones 


B.  R. 


Tcos  6-2  sin# 
c c 

T(l-sin#  ) 

c 


(40a) 


and  for  flat- face  cones 


B.  R. 


Tcos  6-2  sin# 
c c 

T cos  6 

c 


(40b) 
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The  results  given  in  the  remainder  of  this  report  are  presented  as  a 
function  of  the  thickness  ratio.  Through  the  use  of  Eqs.  (40a)  and  (40b)  they 
can  easily  be  interpreted  in  terms  of  the  bluntness  ratio. 

A.  SPHERE-CONES 

The  body  geometries  for  a family  of  10-deg  sphere-cones  with  constant 
base  radius  are  shown  in  Fig.  3a.  For  T ^ 2.1,  the  shapes  are  truncated 
spheres  and  for  T < 2.  1 the  shapes  are  sphere-cones  with  bluntness  ratio 
decreasing  as  T decreases.  As  the  bluntness  ratio  B.  R.  approaches  zero, 
the  thickness  ratio  approaches  the  Sharp  cone  limit 


T , =2  tan#  = 0.  332 

sharp  cone  c 


(41) 


Equation  (41)  follows  directly  from  Eq.  (40a)  for  B.R.  = 0. 

The  relationship  between  the  integrated  heat  transfer  rate  dQ/dt,  the 

freestream  parameters  p , V and  H , the  base  radius  r^,  and  the  heat 

L °o’  oo  o'  B' 

transfer  integrals  Iq  and  Iq  is  shown  in  Fig.  1 and  is  given  by  Eqs.  (10a) 

and  (10b).  If  the  freestream  parameters  and  the  base  radius  are  fixed, 

then  the  relative  heat  transfer  rates  for  different  bodies  can  be  evaluated  by 

knowing  Iq  as  a function  of  thickness  ratio  T. 

The  Iq  as  a function  of  T for  the  family  of  10-deg  cones  in  Fig.  3a  and 
all  other  10-deg  cones  is  shown  in  Fig.  3b.  For  truncated  spheres,  T > 2.  1, 
Iq  is  approximately  constant.  For  sphere-cones  T<  2.  1,  Iq  decreases  as  T 
decreases  due  in  large  part  to  the  decreasing  area  associated  with  the  high 
heat  transfer  spherical  cap  and  the  fact  that  the  inviscid  flow  along  the  cone 
is  overexpanded.  As  T continues  to  decrease,  however,  a minimum  value  of 
Iq  is  reached.  Further  decrease  in  T results  in  an  increase  in  heat  transfer 
because  the  conical  region  now  includes  flow  recompression  following  the 
local  overexpansion. 

A sharp  cone  represents  the  smallest  thickness  ratio  attainable  for  a 
given  cone  half-angle.  The  10-deg  sharp  cone  value  of  Iq  for  isentropic 
expansion  of  the  flow  is  shown  by  an  x on  Fig.  3b. 
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The  blunt  cone  calculated  results  are  seen  to  be  approaching  the  sharp 
cone,  isentropic  expansion  limit  as  T approaches  its  lower,  sharp  cone  limit. 

All  results  presented  in  the  remainder  of  this  report  assume  isentropic 
expansion  of  the  flow  and  are  therefore  applicable  only  when  entropy  swallow- 
ing effects  are  negligible.  This  is  generally  an  excellent  assumption  for 
nosetip  calculations. 

Calculated  values  of  Iq  as  a function  of  T for  laminar  and  turbulent 
flow  over  sphere-cones  are  shown  for  various  cone  half-angles  in  Figs.  4a 
and  4b,  respectively.  The  results  are  generally  similar  to  those  described 
in  Fig.  3.  For  each  cone  half-angle  greater  than  zero,  the  laminar  and  1 

turbulent  curves  of  Iq  versus  T pass  through  a minimum.  The  minimum 
Iq  or  ^ decreases  as  6^  decreases. 

To  illustrate  the  effect  of  cone  half-angle,  bluntness  ratio,  and  thick- 
ness ratio  on  the  heat  transfer  rate,  the  L-.  behavior  illustrated  in  Figs.  4a  < 

^ J 

and  4b  has  been  reproduced  in  Fig.  5a.  The  different  symbols  on  the  curves,  t 

marked  by  numbers,  represent  families  of  cones  to  be  considered  in  the 
following  discussion. 

First  consider  the  four  • (filled  circle)  symbols  numbered  1,  2,  3 and  4. 

Each  point  represents  the  minimum  in  the  curve  of  Iq  versus  T for  a particular  * 

cone  half-angle.  Four  cones,  having  a constant  base  radius  and  with  cone 

half-angles,  thickness  ratios  and  bluntness  ratios  corresponding  to  points 

1 through  4,  are  shown  in  Fig.  5b.  Since  the  base  radius  is  constant,  the 

intercepted  freestream  total  energy  flux  is  constant.  However,  the  fraction  t 1 

of  this  flux  transferred  as  heat  to  the  body  decreases  monotonically  as  the 

cone  half-angle  and  thickness  ratio  decrease,  since  dQ/dt  is  proportional  to 

Iq;  thus 


dQ4  dQ3  dQ  dQ 

~dT>  ~dT>~dT>~dT  <42> 


Iq  as  a Function  of  Thickness 
Ratio  for  Sphere-Cones 


Fig.  5c.  Sphere -Cones  with  Constant  Thickness 
Ratio,  T = 1.  3 


Fig.  5d.  Sphere -Cylinders 
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The  four  o (open  circle)  symbols  in  Fig.  5a  numbered  5,  6,  7 and  8 represent 
four  cones  having  a constant  thickness  ratio  of  1.3.  These  cones  are  illus- 
trated in  Fig.  5c.  The  20-deg  cone,  represented  by  point  number  8,  has  the 
lowest  heat  transfer  rate.  The  5-  and  10-deg  cones,  points  6 and  7,  have 
higher  heat  transfer  rates  because  of  larger  areas  associated  with  the  high 
heat  transfer  nose  cap.  The  30-deg  cone,  point  5,  also  has  a higher  heat 
transfer  rate  because  of  the  high  conical  pressure  on  the  cone  surface.  Thus, 
for  these  cones 


dQ_  dQ  dQ7  dQ 

_ > > > 

dt  dt  dt  dt 


(43) 


The  behavior  of  Iq  versus  T foi;  a sphere-cylinder  can  now  be  under- 
stood by  considering  Fig.  5d  and  the  o (half-filled  circle)  symbols  in 
Fig.  5a,  points  9,  10  and  11.  We  see  in  Fig.  5d  that  body  9 may  be  considered 
a part  of  body  10.  Likewise,  bodies  9 and  10  may  be  considered  parts  of 
body  11.  Therefore,  going  from  body  9 to  body  10  and  then  to  body  11,  the 
total  heat  transfer  must  increase  or 


dQll  dQlO  dQ9 
dt  > dt  > dt 


(44) 


B.  FLAT-FACE  CONES- -COMPARISON  WITH  SPHERE-CONES 

Calculated  values  of  Iq  for  laminar  and  turbulent  flow  over  flat-face 
cones  are  shown  for  various  cone  half-angles  in  Figs.  6a  and  6b,  respectively. 
The  behavior  is  comparable  to  that  for  sphere-cones  as  shown  in  Figs.  4a 
and  4b  with  some  notable  differences.  The  surface  pressure  distributions 
for  a 10-deg  sphere-cone  and  a 10-deg  flat-face  cone  are  compared  in  Fig.  7. 
Since  Iq  depends  upon  the  pressure  distribution  through  Eqs.  (38)  and  (39), 
differences  in  the  pressure  distributions  will  be  reflected  in  the  Iq  versus 
T behavior.  The  major  difference  between  the  sphere -cone  and  the  flat-face 
cone  pressure  distributions  is  in  the  region  of  x/rR  < 1.  0.  The  rapid 
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expansion  of  the  flow  at  the  corner  of  the  flat-face  cone  results  in  a severe 

local  overexpansion  with  immediate  recompression  to  a pressure  level 

greater  than  that  for  the  same  station  on  a sphere-cone.  For  x/r  >1, 

n 

the  sphere-cone  and  flat-face  cone  pressure  distributions  are  similar,  with 
the  flat-face  pressure  remaining  above  the  sphere -cone  value. 

For  the  flat-face  cones,  the  increasing  pressure  with  increasing  x/r^ 
for  x/rn  < 1.0,  is  responsible  for  the  increasing  Iq  with  decreasing  T 
behavior  for  1 < T < 10  as  seen  in  Figs.  6a  and  6b.  For  T < 1,  the  sphere- 
cone  and  flat-face  cone  results  are  qualitatively  the  same. 


The  numerical  values  of  Iq  for  sphere-cones  versus  flat-face  cones 

are  compared  for  laminar  and  turbulent  flow  in  Figs.  8a  and  8b,  respectively. 

In  general,  for  T > 2,  (L_)r1  . - < (I-,)  , and  for 

Q flat-face  cone  Q'sphere-cone 

T < 2’  ^sphere-cone  < ^flat-face  cone*  The  differen«s  are  not  large. 
For  T > 2,  the  maximum  increase  in  Iq  for  a sphere-cone  compared  to  a 

flat-face  cone  is  about  50  percent.  For  T < 2,  the  flat-face  cone  Iq  values 

are  5 to  30  percent  greater  than  the  sphere-cone  values.  From  the  above 

discussion,  the  reasons  for  this  comparative  behavior  are  apparent.  The 

relatively  low  heat  transfer  on  a flat-face  and  the  very  low  local  pressures  on 

the  forward  part  of  a flat-face  cone  account  for  an  Iq  less  than  that  for  a 

sphere-cone  with  the  same  thickness  ratio.  However,  as  x/r  increases, 

which  means  T decreases,  the  higher  flat-face  cone  local  pressure  values 

eventually  result  in  (I-.),.  . . becoming  greater  than  (I_)  . 

Q flat -face  cone  6 6 ' Q'sphere-cone. 

C.  EFFECT  OF  y AND  U 


All  results  presented  thus  far  have  been  for  the  nominal  values  of 
specific  heat  ratio  y and  viscosity  exponent  U)\  i.  e.,  y « 1.2  and  u » 1.0. 

While  these  values  appear  to  be  the  most  appropriate  constant  values  for  the 
selected  « 20  flowfield,  the  question  arises  as  to  the  sensitivity  of  the 
results  to  changes  in  these  nominal  values.  To  illustrate  the  sensitivity, 
calculated  results  are  presented  in  Fig.  9 showing  Iq  for  10-deg  sphere-cones. 
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Fig.  9. 


Sensitivity  of  I^to  Y and  0), 
10-deg  Sphere-Cone 
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for  all  combinations  of  7 = 1.  2 or  1.4  and  0)  = 1.  0 or  0.  7.  Both  increasing  y 
from  its  nominal  value  and  decreasing  (J  from  its  nominal  value  result  in  an 

flj  rm 1 *1 tL 


The  highest  Iq  value  for  given  T occurs  for  y*  1.4 


/ 2 


increase  in  I 
and  = 0.7. 

The  heat  transfer  rate  dQ^/dt  is  proportional  to  ^Iq^kcoJ  as  given 
by  Eqs.  I 

bining  Eqs.  (10a),  (20a)  and  the  results  of  Fig.  9,  we  conclude  that 


['qAo]1 


(dQ  /dt)y=1,4  ^ = 0.  7 
(dQ  /dt)y=l.  2,W=1.  0 


= 1.576 


(45a) 


max 


Similarly 


(dQ  /dt)y-it  4>(j=o>7 

(dQ  /dt)y  = l.  2,W  = 1 


= 1.  161 


(45b) 


max 


I 

> 


I 


| 
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IV.  CALCULATED  RESU LTS-- TRAJECTORY  ENVIRONMENT 


A.  LAMINAR  OR  TURBULENT  FLOWS 

In  this  case,  the  total  heat  transfer  is  given  by  Eqs.  (31a)  and  (36a) 

for  laminar  and  turbulent  flow,  respectively.  From  these  equations,  we  see 

that  for  specified  reentry  conditions  and  vehicle  ballistic  coefficient  the  total 

heat  transfer  Q depends  on  the  body  shape  and  the  body  scale  through  the 

parameters  Iq  and  r^,  respectively.  Since  the  Iq's  in  Eqs.  (31a)  and  (36a) 

are  the  same  as  those  in  Eqs.  (10a)  and  (10b),  results  discussed  above  in 

conjunction  with  Figures  3 through  9 are  also  directly  applicable  to  the 

L t 

present  case.  Thus,  Figs.  4 and  6 supply  the  values  of  Iq  and  Iq  to  be  used 
in  Eqs.  (31a)  and  (36a)  to  evaluate  Q*-1  and  for  sphere-cones  and  flat- 
face  cones,  respectively.  Also  needed  are  the  values  of  the  functions  con- 
taining  the  ballistic  parameters  B and  B , i.e.,  erfvB  and  y(4/5,  B ). 
These  functions  are  shown  in  Fig.  10. 

If  the  calculated  results  given  in  Figs.  4,  6,  and  10  are  used  in  com- 
bination with  Eqs.  (31a)  and  (36a),  the  total  heat  Q[Btu]  transferred  to  a nose 
tip  or  a reentry  vehicle  for  a complete  reentry  trajectory  can  be  easily  cal- 
culated. The  appropriate  units  for  the  parameters  in  Eqs.  (31a)  and  (36a)  are 
CJJ0.979  X 10-8  (lbm/ft-sec-°R)],0[lbm/ft2],  A[ft_1],  rB[ft],  VE[ft/sec]. 

If  the  nose  tip  heat  transfer  is  desired,  then  the  r^  used  in  these 
equations  should  be  the  "base"  radius  of  the  nose  tip,  and  the  thickness 
ratio  Tused  to  evaluate  Iq  from  the  figures  should  be  that  for  the  nose  tip. 
Likewise,  for  evaluating  the  total  heat  transferred  to  the  entire  reentry 
vehicle,  the  r^  and  T should  correspond  to  the  entire  vehicle  values  of  these 
parameters. 

Caution  must  be  applied  in  using  the  results  calculated  in  the  preceding 
manner  if  the  bluntness  ratio  B.  R.  = r^/r^  of  the  nose  tip  or  vehicle  is  less 
than  0.  1,  since  entropy  swallowing  effects  are  then  likely  to  be  important. 


Fig.  10.  Laminar  and  Turbulent  Heat 
Transfer  Functions 


B. 


LAMINAR-TURBULENT  FLOW 


The  above  procedures  allow  the  calculation  of  the  total  heat  transfer  Q 
only  for  those  cases  in  which  the  flow  in  the  boundary  layer  is  all-laminar 
or  all-turbulent  for  the  entire  trajectory.  If,  when  boundary  layer  transition 
occurs,  it  is  assumed  that  the  flow  changes  instantaneously  from  laminar  to 
turbulent  over  the  whole  nose  tip  or  vehicle,  then  the  total  laminar-turbulent 
heat  transfer  Q can  be  evaluated. 

For  the  laminar -turbulent  case  with  boundary  layer  transition  occurring 
at  altitude  h we  have 


The  details  of  performing  these  integrations  are  given  in  Appendix  A.  The 
final  result  is 


QLT  = FLerf 


n/b^+  FT 


y(4/5,BT)  - y(4/5,B^) 


(47) 


where  F^1  and  F^ 
partial  trajectory 


are  defined  by  Eqs. 

ballistic  factors  B^ 
P 


(31b)  and  (36b),  respectively,  and  the 
T 

and  Bp  are  given  by 


L -Xh 
e 

(48a) 

T -Xh 
e 

(48b) 

Thus,  ia  evaluated  in  the  same  manner  as  QL  and  QT  except  multiple 

L 

use  of  Fig.  10  is  required  to  give  the  three  components.  From  either  Q , 
QT,  or  QLr,  the  total  mass  of  an  ablating  material  lost  during  a reentry 
trajectory  is  given  by 


c. 


EVALUATING  LOCAL  HEAT  LOADS 


In  many  instances,  it  is  desirable  to  know  the  total  local  heat  load  at  a 
particular  body  station  in  comparison  to  the  total  body  heat  load.  The  total 
local  heat  load  Jifis  related  to  the  total  body  heat  load  Q by 


= 


1 dQ 


2m*  ds 


(49) 


Expressions  for  dQ/ds  may  easily  be  obtained  from  Eqs.  (31a)  and  (36a)  for 

L T 

laminar  and  turbulent  flow,  respectively.  In  these  equations,  Iq  and  Iq  are 
the  only  variables  which  are  a function  of  s.  Performing  the  indicated  dif- 
ferentiation and  substituting  into  Eq.  (49),  we  obtain  for  laminar  flow 


J * 


1 . 443  x 10 


-7 


mi 


-|l/2 


(2UH3) 


Cp(l^)  P 


Xsin0Er 


1/2 


(3+2o»)/2  / rT7\ 

qLerf(VBL  ) (50a) 


and  for  turbulent  flow 


.X?T  = 1. 691  x 10-7 


C /y+ 1 \ 1 

4/5 

r (y~l  f 

lynv 

^UXZITr)  j 

(2tJ+9) 

[Asin#Er 

1/5 

/J  4/5 


(50b) 
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The  nondimensional  heat  transfer  parameters  q and  q are  given  by 


T p u M r 

L *e  e e 

q = 


WT 


P u P r * ^ 
e e e 

W7^ 


(51a) 


(51b) 
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Comparing  Eqs.  (51a)  and  (51b)  with  Eqs.  (la)  and  (lb),  we  see  that 

T -p 

q and  q are  simply  nondimensional  representations  of  the  original  Vaglio- 
Laurin  laminar  and  turbulent  heat  flux  expressions. 

Li  T 

In  Eqs.  (50a)  and  (50b),  q and  q have  been  normalized  in  such  a way 

L T 

that  q is  unity  at  an  axisymmetric  stagnation  point  and  q is  unity  at  the 

L T 

maximum  turbulent  heating  point  on  a sphere.  The  parameters  q and  q 
have  been  calculated  and  are  shown  in  Figs.  11  and  12  for  sphere-cones  and 
flat-face  cones,  respectively.  Notice  that  the  reference  length  in  Eqs.  (51a) 
and  (51b)  is  the  nose  radius  or  flat-face  step  height,  designated  r^. 

We  are  now  in  a position  in  which  it  is  possible  to  easily  compare  the 
present  results  with  those  of  Brunner.  Recall  that  Brunner's  starting  point 
was  the  flat-plate-reference-enthalpy  heat  transfer  expressions,  whereas  the 
present  starting  point  has  been  the  Vaglio- Laurin  equations.  The  detailed 
algebraic  manipulations  for  this  comparison  are  given  in  Appendix  B.  The 
laminar  total  local  heat  transfer  Jlf1"  given  by  Eq.  (50a)  is  for  an  axisymmetric 
stagnation  point  11.  1 percent  lower  than  that  calculated  from  Brunner's 
expression.  The  turbulent  total  local  heat  transfer  given  by  Eq.  (50b) 

for  the  maximum  turbulent  heating  point  on  a sphere,  is  14.  9 percent  higher 
than  that  calculated  from  Brunner's  equation. 

For  an  ablating  material,  the  total  recession  depth  d at  a given  nose 
tip  or  body  location  can  be  calculated  from 


where  p^  is  the  bulk  density  of  the  ablating  material  and  Q is  the  cold  wall 
effective  heat  of  ablation.  If  Figs.  11  and  12  are  used,  the  d for  any  location 
on  the  families  of  sphere -cones  and  flat-faced  cones  considered  in  this  work 
can  be  easily  calculated. 
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Fig.  12.  Local  Heat  Transfer  Function, 
Fiat-face  Cones 


V.  SUMMARY 


Beginning  with  the  hypersonic  convective  heat  transfer  expressions  for 
laminar  and  turbulent  flow  given  by  Vaglio- Laurin,  equations  for  the  total 
heat  transfer  to  a body  in  a hypersonic  flow  environment  have  been  derived. 
When  these  equations  are  used  together  with  the  graphical  results  presented, 
total  heat  transfer  to  selected  body  geometries  can  be  calculated  directly. 

Both  the  case  of  a constant,  ground  test  type  environment  and  a realistic 
reentry  vehicle  trajectory  environment  have  been  considered. 

The  total  heat  transfer  is  given  as  a function  of  freestream  or  reentry 
parameters,  body  scale  and  heat  transfer  integrals  which  are  functions  of  the 
body  geometry.  The  best  available  inviscid  flow  field  calculation  methods  for 
surface  pressure  distributions  were  used  to  calculate  the  laminar  and  turbu- 
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lent  heat  transfer  integrals  Iq  and  Iq  for  0-,  5-,  10-,  20-,  and  30-deg  sphere- 
cones  and  flat-face  cones.  These  results  are  presented  graphically  as  a 
function  of  the  fineness  ratio,  i.  e.,  the  diameter  to  length  ratio.  From  these 
figures,  the  variation  of  total  heat  transfer  with  nose  bluntness  ratio,  fine- 
ness ratio,  and  cone  half-angle  has  been  examined. 

For  the  trajectory  environment,  in  addition  to  total  laminar  or  turbulent 
heat  transfer  prediction  methods,  a method  for  evaluating  the  total  heat 
transfer  for  laminar  flow  to  boundary  layer  transition  altitude  and  all  turbu- 
lent flow  to  impact  is  also  given.  When  the  total  heat  load  to  the  body  or  nose 
tip  is  known,  the  heat  shield  weight  or  transpiration  nose  tip  coolant  require- 
ments can  be  directly  calculated. 

Finally,  graphical  results  giving  the  distribution  functions  for  local 
laminar  and  turbulent  heat  transfer  over  the  sphere -cone  and  flat-face  cone 
geometries  considered  in  this  work  are  given.  These  results  allow  the  total 
local  heat  load  at  selected  axial  stations,  and  therefore  the  local  depth  of 
heat  shield  ablation,  to  be  calculated. 
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DERIVATION  OF  QLT  RELATIONSHIP 
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For  a laminar-turbulent  trajectory  case,  Q is  evaluated  by  per- 
forming the  integrations  indicated  by  Eq.  (46).  The  integrand  dQ^/dy  is 
given  by  Eq.  (26a)  and  dQ  /dy  is  given  by  Eq.  (32a).  These  integrands 
have  already  been  integrated  between  the  limits  of  zero  and  infinity  with  the 
final  results  given  by  Eqs.  (31a)  and  (36a).  We  now  consider  the  partial 
integrations . 

First,  for  the  laminar  portion  of  the  trajectory 


h 


rL  -Ay/ 2 -BLe'Xy, 
C,  e ’ e dy 


(A  - 1) 


When  the  substitution  t = \T&L  e Xy^2  is  made,  Eq.  (A-l)  may  be  rewritten 

7*1 


as 


-c3l^ 


Qp  " [BL]1/2A 


e-t2dt 


(A-2) 


Equation  (A-2)  is  the  same  form  as  Eq.  (31c)  except  that  the  upper  limit  of 
integration  is  different.  Thus,  we  define  B^1  by 


BL  = BLe'Xh 
P 


(A- 3) 


The  integral  in  Eq.  (A-2)  is  thus  represented  by  the  error  function  of  argu- 
ment \/ B ^ and  the  first  term  in  Eq.  (47)  is  obtained  after  appropriate  sub- 
stitutions. 


For  the  turbulent  portion  of  the  trajectory 


T -(4A/5)y  -rB  e 
e e 


T -Ay 


(A-4) 


When  the  substitution  t = is  made,  Eq.  (A-4)  becomes 


T [BT]4/5 

iT  = ~5C3  r 

P ~ 4\[BTj4/5  J[BT]4/5e-(4M5)h 


5/4 

e dt 


(A-5a) 


4A[BT]4/5 


iRT14/5 
B 1 .5/4 

e'  dt 


■rT-i4/5  - (4A/5)h 

B J e .5/4 

e“  dt 


(A-5b) 


The  first  term  of  Eq.  (A-5b)  is  the  same  as  Eq.  (37).  The  second  term  dif- 

T 

fers  only  in  the  upper  limit  of  integration.  Thus,  we  define  B by 


BT  = BTe"*h 
P 


(A -6) 


Then  with  the  definition 


■(4/5,  x)  = I 


(A -7) 


the  last  two  terms  in  Eq.  (47)  follow  directly  from  Eq.  (A- 5b)  after  appro- 
priate  substitution. 
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APPENDIX  B 


COMPARISON  WITH  BRUNNER'S  WORK 


The  major  thrust  of  the  present  work  has  been  to  obtain  Eqs.  (31a) 

Li  T 

and  (36a),  which  give  the  total  heat  transfer  values  Q and  Q . These  equa- 
tions were  obtained  by  performing  successive  body  surface  and  time  (altitude) 
integrations  of  the  appropriate  local  heat  transfer  rate  and  altitude -dependent 
expressions  for  heat  transfer.  Brunner's  work  involved  only  the  time  (alti- 
tude) integrations  of  the  local  heat  transfer  rate  to  obtain  the  total  trajec- 
tory integrated  heat  transfer  at  a point  on  the  body. 

As  indicated  by  Eq.  (49),  d and  Q are  related  by 


„„  dQ/ds 
^ = 2mr 


(B-l) 


Since  the  starting  point  of  the  present  work  is  the  Vaglio- Laurin  con- 
vective heat  transfer  relationships,  and  Brunner's  starting  point  was  the  flat- 
plate  reference  enthalpy  heat  transfer  expressions,  it  is  of  interest  to  com- 
pare the  values  of  dobtained  from  the  two  methods. 

For  convenience,  we  define  the  following: 

- devaluated  using  Brunner's  equations 

•*2  -M  evaluated  using  the  present  equations 


The  and  df  are  given  by  Eqs.  (50a)  and  (50b),  respectively.  If  the 
present  notation  is  used,  Brunner's  expressions  are  given  by 

f 1 1/2 


4. 155  x 10 


-10 


nfi 


1.  lAsin0ErN 


2.  392  x 10 


-10 


A 


,4/5 


Xsin0ErN 


T7T 


V2’2  erf  (/TTb) 

(B-2a) 

2,48  y(4/5,  1.24B) 

(B-2b) 
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A direct  numerical  comparison  of  the  two  methods  is  now  possible. 

Equations  (B-2a)  and  (B-2b)  are  for  an  axisymmetric  stagnation  point  and  the 

maximum  turbulent  heating  point  on  a sphere,  respectively.  Thus,  we  set 
L T 

q and  q in  Eqs.  (50a)  and  (50b)  to  unity.  Note  also  that  the  exponents  on 
p,  x,  sin^E  and  rN  are  the  same  in  Eqs.  (B-2a)  and  (50a)  for  laminar  flow 
and  in  Eqs.  (B-2b)  and  (50b)  for  turbulent  flow.  Forming  the  ratio 
and  substituting  from  these  equations,  we  obtain  for  a laminar  stagnation 
point  on  a sphere 


1.61 
1.  811 


0.  889 


(B-3a) 


and  for  the  turbulent  maximum  heating  point  on  a sphere 

T 

**2  . 12.75  _ , 

~7t  TT7T  " 1’ 149  (B”3b) 

**1 

For  these  numerical  evaluations,  the  following  input  quantities  were  assumed: 

y = 1.2 

(J  =1.0 


